Ultracold lattice gases with periodically modulated interactions 
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We show that a time-dependent magnetic field inducing a periodically modulated scattering length 
may lead to interesting novel scenarios for cold gases in optical lattices, characterized by a nonlinear 
hopping depending on the number difference at neighboring sites. We discuss the rich physics 
introduced by this hopping, including pair superfluidity, exactly defect-free Mott-insulator states for 
finite hopping, and pure holon- and doublon superfluids. We also address experimental detection, 
showing that the introduced non-linear hopping may lead in harmonically trapped gases to abrupt 
drops in the density profile marking the interface between different superfluid regions. 

PACS numbers: 37.10.Jk, 67.85.Hj, 73.43.Nq 



Ultracold atoms in optical lattices formed by laser 
beams provide an excellent environment for studying 
lattice models of general relevance in condensed-matter 
physics, and in particular, variations of the celebrated 
Hubbard model [l|, Q ■ Cold lattice gases allow for an 
unprecedented degree of control of various experimental 
parameters, even in real time. In particular, intcrpar- 
ticle interactions can be changed by means of Feshbach 
resonances Moreover, recent milestone achievements 
allow for site-resolved detection, permitting the study of 
in-situ densities 0, HI, and more involved measurements, 
as that of non-local parity order Q. 

The modulation of the lattice parameters in real time 
opens interesting possibilities of control and quantum en- 
gineering. In particular, a periodic lattice modulation 
translates by means of Floquet theorem 0] into a mod- 
ified hopping constant Q, which may even reverse its 
sign as shown in experiments 0, [l^l ■ This technique has 
been employed to drive the Mott-insulator (MI) to super- 
fluid (SF) transition ll|, and to simulate frustrated clas- 
sical magnetism [1^ . Recent experiments have explored 
as well the fascinating perspectives offered b y p eriodically 
driven lattices in strongly correlated gases [lil . [l^ . 

The effective Hubbard-like models describing these ul- 
tracold lattice gases are typically characterized by a hop- 
ping rate which is independent of the number of particles 
at the sites. This is, however, not necessarily the case. 
Multiban d p hysics may lead to occupation-dependent 
hopping 15|, ll6|. In addition, long-range dipole-dipole 
interactions may lead to number-dependent hoppings as 
well, for sufficiently large dipole strengths A major 
consequence of non-linear hopping is the possibility to 
observe pair superfluidity (PSF) [l3, [ill- 

In this Letter, we consider a cold lattice gas in the pres- 
ence of a periodically modulated magnetic field. In the 
vicinity of a Feshbach resonance, this field induces mod- 
ulated interparticle interactions. Interestingly, Ref. [l^ 
has shown that such periodic modulations of the inter- 
action potential may lead to a many-body coherent de- 
struction of tunneling in two- mode BECs. As shown 
below, the generalization of this effect to lattice gases. 



leads under proper conditions to an effective Hubbard- 
like model with a non-linear hopping which, in contrast 
to other proposals mentioned above, depends on the dif- 
ference of occupations at neighboring sites, and retains its 
non-linear character even for weak lattices. We discuss 
the rich physics introduced by this hopping, including 
PSF phases, exactly defect-free MI states for finite hop- 
ping, and pure holon- and doublon superfluids. We also 
address experimental detection, showing that the stud- 
ied non-linear hopping may lead to abrupt drops in the 
density profile of harmonically trapped gases. 

We consider bosons in a lattice in the presence of a 
periodically modulated magnetic field B{t) = B{t + T) 
(with period T = 2Tr/u!) chosen close to a Feshbach res- 
onance, where the s-wave scattering length acquires the 

form a{t) = a^g {l + B(f)^B^ ^ = ao + Y.iyo^i'^'^'^i^^^)- 
Here A_B and determine the width and position of 
the resonance, respectively, and a^g is the background 
scattering length [3|. Assuming that the gap between 
the first two lattice bands is much larger than any other 
energy scale in the problem, we consider only the low- 
est band and describe the system by a Bose-Hubbard 
model (BHM) 0,0: 



U{t) 



^ rii (n, - 1) - ^ n,, (1) 



where hi {b\) is the bosonic annihilation (creation) op- 
erator at site i, fii = h\hi, /x is the chemical poten- 
tial, J > is the hopping rate and (..) denotes nearest 
neighbors. Interactions are characterized by a coupling 
Uit) = Uo + T,i>o Ui cos(/wi) = Uo + U{t), with Uq > 
and Ui = ^^^^ J d^r\w{r)\^. Here w{r) is the lowest 
Wannier function and M is the atomic mass. 

We apply a similar analysis as the one used for shaken 
lattices [8|. We specify a Floquet basis 



\{nj},', 



, (2) 



where m defines the Floquet sectors and is the 

Fock basis, characterized by the atom number at each 
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FIG. 1: (Color online) Phase diagram for — 4 using 
GA (left) and DMRG (right). Solid curves define the MI 
lobes, whereas dashed curves are the SF-PSF boundaries. 

site. In Eq. we defined V(t) = ftj{t')dt'/h. 

We introduce the time-averaged scalar product 

(({n^},m'|...|{rij},m)) = m'| . . . m), 
and establish the matrix elements: 

{{{n'^},m'\[H{t)~ihdt]\{n,},m)) 

= Sm.rn' ({"-j } I I {"j }> 

hi ~ hj)bj\{nj}), (3) 

with Hm = mhuj + hj {•fij — 1) — ^fij and 

F,n{x) = ^ /J'dte^""*e'^(*)^. If fvjj > J,Uq, we may 
restrict to a single Floquet sector m = 0, resulting in an 
effective time-independent Hamiltonian of the form 

-ffeff = -jy^^hlFgifii - nj)bj 

+ T ^ ^^^^ i-^Yl 

i i 

Hence, interactions with a periodic modulation result in 
a non-linear hopping term, which depends on the atom 
number difference between neighboring sites. Note that 
this non-linear character remains relevant for any value 
of the bare hopping J. In the following we discuss the 
specific case U{t) ~ Ui cos ujt. In this case, Fq{x) = 
Jo (^2;), with Jo the Bessel function and ft = Ui/huj, 
generalizing the result of Ref. [l^ for two-well BECs. 

Insight on Eq. @ is gained by means of a Gutzwiller 
Ansatz (GA) for the ground state \G) = 

rij Z^n where /„(j) are variational param- 
eters (X]nl/™(j)P ^ l)j determined by minimizing 
{G\Hos\G). Results by choosing homogeneous real 
= fn dH are shown in Fig. [T] (left), where we depict 
the mean- field phase diagram for fi = 4 (Jo(^^) — —0.4) 
as a function of fi/Uo and zJ/Uq, with z the coordi- 
nation number. As usual, MI phases are characterized 



by integer {hi), and vanishing single- particle- and pair- 
condensation parameters, {bi) and {bf). In addition to 
the usual SF phase, we find a PSF phase characterized 
by > \{bi)\'^- Pair superfiuidity is especially pro- 

nounced in the vicinity of integer (hi). Our GA results 
show that PSF phases only occur if Jo(^^) < 0. This may 
be understood by considering integer filling (hi) = n, and 
restricting the variational space to /„±i = ^^^e"''* and 
/„ = cos 77. For Jo(f^) < 0, energy minimization gives 
Alp = — ip- = TT, while Aip = for Jo{il) > 0. As a 
result, for dominant hopping {2Jz/Uq ^ 1), PSF needs 
A^n{n 4- 1) > (^ + sign( Jo(^^))\/'^ + 1)^, which can 
only be fulfilled if Jo(f}) < 0. 

To complement the mean-field GA results, we have 
also employed numerically exact methods in ID. In 
particular, we used the density-matrix renormalization 
group (DMRG) |22| with up to 40 sites and keeping 200 
states, and a related method, the infinite time-evolving 
block decimation (iTEBD) [23[ using a Schmidt dimen- 
sion of 200. We have monitored the behavior of single- 
particle and pair correlations, Gi{i,j) = {b\hj) and 
G2{i,j) = {{^ly^b'j), respectively. Both decay exponen- 
tially in the MI. PSF phases are characterized by a dom- 
inant G2, i.e., a slower power-law decay than Gi. The 
opposite characterizes the SF phase. Figure [T] (right) 
shows the ID phase diagram for J7 = 4, which closely 
resembles the GA one. Similar to the GA, we observe a 
PSF phase, which for integer (h) approaches all the way 
to the tip of the MI lobes. On the other hand, away from 
the lobe tips we observe a direct MI-SF transition. Our 
ID results also confirm the absence of PSF for Jo(^) > 0. 

The case Jo{fl) = is particularly interesting, since 
for neighboring sites i and j with equal number of parti- 
cles, the process \n)i\n)j — >■ jnib l)i|n=F 1)^ is forbidden. 
However, the hopping |n ± l)i\n)j — > |n)i|n± 1)^ is still 
characterized by the usual rate J. This difference has a 
remarkable impact for both the MI- and the SF phases. 

For J = 0; the ground state of Eq. ^ is. as for the 
standard BHM (fi = 0), a defect- free MI {g)^- \n)j for 
n - 1 < n/Uo < n For and J > 0, this 

state is not an eigenstate of Eq. ^ , and quantum fluctu- 
ations induce a finite particle-hole population in the MI 
with an associated non-local parity order Q. Interest- 
ingly, the defect-free state remains an eigenstate of 
for Jb(^^) = 0. As a result, the whole MI lobe is charac- 
terized by the absence of particle-hole defects. Although 
this is typically an artifact in the mean-field GA, in this 
case it is an exact result, and indeed our ID iTEBD 
results show a vanishing variance (An)^ = (n^) — (h)'^ 
within the Mott region (see Fig. [3]). 

Conversely, particles or holes (|n ± 1)) on top of the 
state \n)j acquire also remarkable properties. For 
any Q, extra particles and holes move with a hopping 
rate (n -I- 1) J and nJ, respectively. For $7 = defects 
are unstable, being created and destroyed by processes 
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J/U^ 

FIG. 2: (Color online) iTEBD results for the holon (circles) 
and doublon (squares) populations as a function of J/f/o for a 
ID system with ^/Uo = 0.3 and = 2.405. Note the absence 
of defects in the MI {J/Uo < 0.15), and the appearance of the 
holon (H-SF) and doublon SF (D-SF). 

\n)i\n)j o |n ± l)i\n =F Since these processes are 
forbidden for Joi^) = 0, defects remain stable [2^. Ne- 
glecting occupations other than n and n ± 1, the defects 
are described by an effective Hamiltonian Hh+Hp, where 

Hh = -JnJ2hlhj + {ji-Uo{n-l))Y,hlh,, (5) 

<i,j> i 

Hp = ~J{n + l)Y,plPJ + iUon-^l)Y,plp^^ (6) 

<ij> i 

characterize, respectively, the physics of holes and parti- 
cles, with the hard-core assumption plpi + hlhi = or 1, 
with hi {pi) the bosonic operators for extra holes (parti- 
cles) at site i. In Eqns. ([5]) and ^ we have set the energy 
of the defect-free MI state Eui = 0. Thus the system be- 
haves as a two-component hard-core lattice Bose gas. For 
higher dimensions, a dilute gas of extra holes (holon gas) 
may be considered as a basically free (supcrfluid) Bose 
gas, with a dispersion £'/i(q) = — L/o(n — 1) + ^6°, where 
e° = — 2JV ■ ^ ^ cosiqid) for a 3D cubic lattice and d 
is the lattice spacing. On the other hand, the dilute gas 
of extra particles ("doublon" gas [2^) has a dispersion 

£^p(q) = [/o'^-Ai+ (n + l)e°. 

At zero temperature, the defect gas condenses for /i < 
/ic = U^in — 1/2) — Jz at the bottom of the holon band, 
i?;i(0), acquiring a pure holon character. On the other 
hand, for [i > [i^ the system condenses at -E'p(O) into 
a pure doublon gas. Hence, remarkably, we expect an 
abrupt jump of (n) (i.e. a diverging compressibility) at 
the line /i = /ic, which coincides with the line of integer 
(71) = n. Figure [3] depicts our GA results for the density 
as a function of [ijV^ and J /Uq, which, as expected from 
the previous discussion, presents an abrupt jump between 
a holon and a doublon superfluid. 

In ID, the defects behave, due to the hard-core con- 
straint, rather as a two-component Tonks gas, but a simi- 
lar two-band reasoning as above applies, and we may also 
expect the existence of pure holon and doublon supcrflu- 
ids. Figure [2] shows our iTEBD results in the vicinity of 




FIG. 3: (Color online) Homogeneous GA results for (n) as a 
function of J/C/q and f^/Uo for Q, = 2.405. Red curves denote 
the boundary of the MI and the line of integer filling 1. Note 
the abrupt jump in the density at that line, indicating the 
transition between the holon SF and doublon SF regimes. 

{n) = 1 for the holon (doublon) populations (Pq) ((P2)) 
, with Pn = Yin'^ni''^ ~ '^')/("- ~ "■')• addition to the 
MI phase characterized by (Pq) = {P2) = 0, wc observe 
a holon-SF ((P2) = 0) and an abrupt jump to a doublon- 
SF ((Po) = 0). Note that pure doublon- or holon-SF 
exclude PSF. 

At constant ^ the system undergoes a MI - dou- 
blon (holon) SF transition at a critical tunneling J dp) 
for which Ep(^h){^) = Emi- In ID, our iTEBD 
results show that this transition retains a commen- 
surate/incommensurate nature as in the usual ID 
BHM [ill, characterized by (Po,2) ^ \/J - Jc- This 
growth is illustrated for the Mott - holon-SF transi- 
tion in Fig. [2] On the other hand, at constant inte- 
ger (fi), there is no ID MI-SF transition at finite hop- 
ping J. This result, expected from the theory of two- 
component Tonks gases [28| , is due to the absence of pro- 
cesses \n)i\n)j o |n ± l)i|n=F 1) j which precludes that 
doublons and holons can swap their positions through 
second-order super-exchanges. As a result, if holons and 
doublons coexist, (which only happens at the singular 
integer filling line) superfluidity is absent. Our DMRG 
results for (n) = 1 confirm indeed that Gi^2 decay expo- 
nentially for any finite J. Interestingly, there is however 
a clear transition between a defect-free insulator and an 
insulator with a finite density of holon-doublon pairs. 

For a finite but small Jq{VI), the SF regions retain to a 
large extent their holon/doublon character, although the 
concentration of doublons/holons in the holon/doublon 
SF increases for growing Jq{^) and J. The coexistence 
region for holons and doublons is hence not any more 
singular, although it remains characterized by a large 
compressibility for small Jq{U,). For Jq{VI) < this co- 
existence region becomes the PSF phase discussed above. 
Away from the Mott-tip a direct MI-SF transition is ob- 
served, as discussed above, since at the MI boundary 
holons and doublons do not coexist. 

Let us finally discuss some experimental questions. 
Optimal experimental conditions for periodically mod- 
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— <n(x)> (LDA) <> <n(x)> (trap) 

- <P2(x)> (LDA) □ <P,(x)> (trap) 
■ <Po(x)> (LDA) o <Po(x)> (trap) 
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FIG. 4: (Color online) GA results for the site densities {fij) 
and the on-site holon and doublon populations for a 2D lattice 
with a harmonic confinement V{jx,jy) = Vo{jx + jy)j with 
Vo/J — 0.0075, interaction Uo/J = 5.33, a central chemical 
potential /lo/J = 3, and Q = 2.45 {Jo{0.) ~ 0). Note the 
central doublon-SF region, surrounded by a holon-SF ring, 
and the abrupt density drop separating both regimes. Lines 
indicate local-density approximation (LDA) results. 



ments at selected holding times with V{t) = provide, 
however, a direct conversion between both basis. Inter- 
estingly, time averaged measurements over times much 
larger than the period T just provide the contribution of 
|n, n ± 1) — >■ |n ± 1, n) processes. 

In summary, periodically modulated interactions lead 
to a rich physics for cold gases in optical lattices, charac- 
terized by a nonlinear hopping depending on the number 
difference at neighboring sites. This hopping can lead 
to pair supcrfluid phases, and also to defect-free Mott 
states, and holon- and doublon superfluids, which may 
be revealed by abrupt jumps of the in-situ densities in 
harmonically trapped lattice gases. 

We acknowledge financial support by the Cluster of 
Excellence QUEST. 



ulated interactions are provided by ^^Rb, which has a 
particularly large abg ~ —AOOas (with gb the Bohr ra- 
dius), and a broad Fcshbach resonance at Br ~ 155. 2G, 
with a width AB = -11.6G The desired form 

a(t) ~ flQ + fli cos{ujt) can be achieved for a magnetic 
field dependence B{t)/G ~ 167.56 -I- 5.58 cos(a;t), with 
flo ~ 20aB and ai w 200aB. We consider a lattice spac- 
ing d ~ 0.5/im, and potential depth Vh = sEn, where 
Eji = H'^tt'^ /2'mcP is the recoil energy. For s « 17 
(J < Uq), the value n = 2.4 (Jo(f^) ^ 0) is obtained 
for w ~ 27r X 900 Hz > [/o/^ = Stt x 217 Hz, ensuring 
that only one Floquet manifold is relevant. 

In order to address the question of detection, we have 
to consider the transformation ([2]) between the Floquet 
\{nj},m) and the Fock |{7T.j}) basis. The densities {fii) 
are equivalent in both, therefore the large compressibil- 
ity regions characteristic of |j7o(i^)| — may be revealed 
in in-situ experiments with an additional harmonic con- 
finement. This is illustrated in Fig. 21 where we show 
inhomogeneous GA results for a harmonic trap in 2D. 
As expected from local-density approximation, we ob- 
serve an abrupt jump in the (n,) when the local chemical 
potential crosses its critical value. Note the wedding- 
cake-like form of the profile, characterized by a central 
doublon-SF, surrounded by a holon-SF ring. 

The interpretation of measurements of other observ- 
ables, as e.g., the momentum distribution in time-of- 
flight (TOF) measurements, is more involved, since 
the basis conversion is non-linear, {{n'j}\blbj\{nj}) ^ 
e-iVWi^-"'^+^){{n'^}^m\blbj\{nj},m). However, for the 
holon- and doublon-SF phases the TOF measurement is 
almost time-independent for small |j7o(fl)|. Indeed this 
weak dependence is in itself a proof of the strong holon 
or doublon character of the SF. For large |Jo(^^)| the 
nonlinear conversion is an issue, and in general mea- 
surement results are periodic. Stroboscopic measure- 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885-964 (2008). 

[2] M. Lewenstein et ai, Adv. in Phys. 56, 243 (2007). 

[3] C. Chin et al. Rev. Mod. Phys. 82, 1225-1286 (2010). 

[4] W. S. Bakr et al, Science 329, 547 (2010). 

[5] J. F. Sherson et ai, Nature 467, 68 (2010). 

[6] M. Endres et ai, Science 334, 200 (2011). 

[7] M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

[8] A. Eckardt, C. Weiss, and M. Hohhaus, Phys. Rev. Lett. 
95, 260404 (2005). 

[9] H. Lignier, et ai, Phys. Rev. Lett. 99, 220403 (2007). 
[10] E. Kierig et al, Phys. Rev. Lett. 100, 190405 (2008). 
[11] A. Zenesini et ai, Phys. Rev. Lett. 102, 100403 (2009). 
[12] J. Struck et ai, Science 333, 996 (2011). 
[13] Y.-A. Chen et ai, Phys. Rev. Lett. 107, 210405 (2011). 
[14] R. Ma et al, Phys. Rev. Lett. 107, 095301 (2011). 
[15] O. Dutta, et ai, New J. Phys. 13, 023019 (2011). 
[16] U. Bissbort, F. Deuretzbacher, W. Hofstetter, preprint, 

[arXiv:1108.6047]. 
[17] T. Sowihski, et al, Phys. Rev. Lett. 108, 115301 (2012). 
[18] K. P. Schmidt, et al, Phys. Rev. B 74, 174508 (2006). 
[19] J. Gong, L. Morales- Molina and P. Hanggi, Phys. Rev. 

Lett. 103, 133002 (2009). 
[20] D. S. Rokhsar and B. G. Kothar, Phys. Rev. B 44, 10328 
(1991). 

[21] Allowing different f„ on the two sublattices of a bipartite 
lattice did not lead to variational improvement. 

[22] U. SchoUwock, Rev. Mod. Phys. 77, 259 (2005). 

[23] G. Vidal, Phys. Rev. Lett. 98, 070201 (2007). 

[24] M. P. A. Fisher, et al, Phys. Rev. B 40, 546-570 (1989). 

[25] Processes such as \n + 2)i\n— l)j •<-> \n + l)i\n)j may still 
occur, but for (n) in the vicinity of n they are much less 
probable, especially at low hopping and low n. 

[26] The term "doublon" is used here as a general concept 
for a dilute gas of extra particles on top of a MI at any 
integer (n). 

[27] T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford University Press (2003). 

[28] S. Chen, J. Cao and S.-J. Gu, Europhys. Lett. 85, 60004 
(2009). 

[29] J. L. Roberts, et al, Phys. Rev. Lett. 81, 5109 (1998). 



